Critical fluctuations in spatial complex networks 
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An anomalous mean-field solution is known to capture the non trivial phase diagram of the Ising 
model in annealed complex networks. Nevertheless the critical fiuctuations in random complex 
networks remain mean-field. Here we show that a break-down of this scenario can be obtained when 
complex networks are embedded in geometrical spaces. Through the analysis of the Ising model on 
annealed spatial networks, we reveal in particular the spectral properties of networks responsible 
for critical fiuctuations and we generalize the Ginsburg criterion to complex topologies. 
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Large attention has been recently addressed to the 
effects that different topological properties may induce 
on the behavior of equilibrium and non-equilibrium pro- 
cesses defined on networks and to the possible implica- 
tions for the study of several social, biological and tech- 
nological networks [l|, Heterogeneous degree distri- 
butions, small world and spectral properties, in partic- 
ular, have been recognized as responsible of novel types 
of phase transitions and universality classes (H-Q- For 
instance, scale-free networks present a complex critical 
behavior for the Ising model, percolation and spreading 
processes, that explicitly depends on the exponent of the 
power-law in the degree distributions [ll-Q . On the other 
hand, the existence of non trivial spectral properties is 
crucial for the stability of synchronization processes and 
0{n) models [i]. 

Despite the large interest in the subject, much smaller 
attention has been devoted to critical phenomena on 
complex networks embedded in a metric space 04^, 
though some important problems related to navigabil- 
ity, efficiency and search optimization in spatial networks 
have been already discussed in the literature [iMl- In 
fact, spatial embedding is a very relevant aspect of in- 
frastructure and technological networks, including air- 
ports networks, the Internet, and power-grid networks. 
Moreover, a pivotal role in shaping the topology of social 
networks is played by hidden metric structures in some 
underlying abstract space, such as that of the social dis- 
tance between individuals [1, 0] . 

The aim of this Letter is to investigate the role of spa- 
tial embedding in relation with the critical behavior of 
phase transitions in complex networks. It is well known 
that in regular lattices, space dimensionality governs the 
critical behavior of equilibrium and non-equilibrium sys- 
tems. In particular, below the upper critical dimension, 
critical fiuctuations that are not captured by the mean- 
field approach set in. Similarly, for complex networks 
embedded in a low dimensional space we can expect that, 
as the link probability becomes short ranged, the effect 
of the underlying space might change the critical behav- 
ior leading to a break-down of the validity of (hetero- 



geneous) mean-field arguments. This should be relevant 
to understand real phenomena in spatial networks, such 
as the spreading of viruses 0, the emergence of con- 
gested phases in the packet-based traffic on technological 
networks and cascading failure phenomena in powe 
r-grid networks 

As a prototypical example of the complex behavior in- 
duced by spatial embedding, in this Letter we consider 
the Ising model on annealed scale-free networks. On 
a scale-free network with a degree distribution P{k) ^ 
^-7SF ^ ^j^g critical temperature of the Ising model di- 
verges for "fsF < 3. The critical exponents, computed by 
means of the annealed network approximation [l6[ or by 
assuming a quenched randomness 17, T^, deviate from 



the mean-field ones as long as "fsF < 5, with the excep- 
tion of 7, 7' describing the divergence of the magnetic 
susceptibility x close to the critical temperature Tc, ( 
X ~ |r — Tc\~^''^ ). In fact 7,7' remain always fixed to 
their mean field value 7 = 7' = !. For these reasons 
we refer to the critical behavior of random scale-free net- 
works as the heterogeneous mean-field solution. We de- 
rive here a Ginsburg criterion [l9[ for spatial complex 
networks determining the condition under which critical 
fluctuations become larger than the ones predicted within 
a mean-field approach. In particular, we will show that 
the critical behavior is always mean-field, whenever the 
matrix p = {pij}i,j=i^...,N, fixing the probabilities of ex- 
istence of each link {i,j) has a finite spectral gap A be- 
tween the maximal eigenvalue A and the second maximal 
one A2 . On the contrary, when the spectral gap A in 
the thermodynamic limit, the critical behavior depends 
on the behavior of the tail of the spectrum of p. We will 
demonstrate by theoretical and numerical results that 
the behavior of such tail is well captured by an exponent 
6s, related to the effective dimension deff of the network 
through the relation Ss = {deff — 2)/2. We find that for 
(5s < 1 the critical fluctuations be come dominant and 
close enough to the critical temperature the mean-field 
theory is not sufficient to correctly characterize the criti- 
cal exponents, possibly calling for renormalization group 
calculations. 
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Networks with spatial embedding - We consider net- 
works of TV nodes embedded in a d-dimensional euclidean 
metric space, each node i = I, . . . , N having position r,. 
The minimal hypothesis (20| that can be made on random 
networks with heterogeneous degrees and spatial embed- 
ding is that links are drawn with probability pij 
given by 



Pij 



9jJ{ri,rj) 



I + 0,6, Jin, f,) 



(1) 



where we assumed that [maxi(0)]^[maxf7^' J(r, f^)] <^ 1 
and that the matrix J(fi,rj) only depends on the dis- 
tance between the nodes, i.e. J{ri,fj) — J{\ri — rj\). 
In this ensemble the degree ki of a node i is a Pois- 
son random variable with expected degree ki fixed by 
means of the hidden variables 9i and given by the relation 
ki = ^jPij- Therefore, given a set of expected degrees 

{ki}, we can evaluate the {9i} variables by solving the 
equations ki = J2j Pij ■ Networks with homogeneous de- 
grees are generated by fixing 9i ~ 6 \/i, that corresponds 
to the Manna-Sen model of spatial networks . Another 
special choice is that of space-independent couplings 

Jy = J Vi, j = l,...,iV, that gives 9, k,/^Jj(k)N 

where (fc) = ki/N. In this case, our formalism easily 
recovers known results for both the percolation thresh- 
old and the critical temperature of the Ising model on 
complex networks without spatial embedding. 

Ising model on annealed spatial networks and the Gins- 
burg criterion - We consider a system of binary spin vari- 
ables Si = ±1, for i = 1, . . . ,N, defined on the nodes of a 
given annealed network with spatial embedding and link 
probability given by the matrix p. The partition function 

[l^ for this problem is given by 



{Si} 



with 



H{s,) = -IY1 s^^iJ^j^jSj - 5Z 



(2) 



(3) 



In order to derive a Ginsburg criterion for this statistical 
mechanics problem, we generalized the classical approach 
by means of stationary phase approximation [l9[. Con- 
sidering only the first-order terms in the expansion leads 
to mean-field results. Thus the validity of the mean-field 
solution can be checked by evaluating the higher order 
corrections at the critical point. Critical fiuctuations that 
are neglected by the mean-field set in when the second 
order corrections diverge dominating the behavior of the 
susceptibility at the criticality. In the stationary phase 
approximation, the magnetization of the system is given 
by the to^'s satisfying the self consistent equations 



tanh[/3(if, + ^0,Jy%7n^")] 



(4) 



At the second order of the stationary phase approxima- 
tion [l9| . performing a Legendre transformation we can 
evaluate the free energy T{{mi}) as 



r {{mi}) ^ - \y^^mi9iJ,j9 



^ Y^[{1 - m,) hi(l - m,) + (1 + m,) ln(l + m,)] 



1 

2lp 



lndet[% - l3J,j9,9j{l - m^)] 



(5) 



where the external field Hi — dr{{mi})/dmi and we 
have introduced the parameter z in order to keep track of 
the different orders in the expansion. The susceptibility 



matrix is defined as 



We compute it in the 



paramagnetic phase, where nii = 0, and then we perform 
the projection along the eigenvector associated to the 
eigenvalue A of the connectivity matrix, obtaining 



< 1 1 V- r 



/3pr\d^'f^ (6) 



where 1 is the identity matrix. The instability of the 
paramagnetic phase is now determined in terms of the 
largest eigenvalue A of the matrix pij through the con- 
dition Xa^(-^c) = 0. If we express the susceptibility in 
terms of the spectral density p{X) of the matrix p as 



XA'(T) = -A + r+i J dXpiX)- 



A 



T 



to leading order in 1 /z the critical temperature Tc is g: 

by 



Te = A - - 
z 



dXp{X) 



1-i 



(7) 



;iven 



(8) 



Using dH), we can express the susceptibility, Eq. ([7]), for 

T ^Tr as 



X 



T-Tr. 



1 



1 



dX- 



p{X){Xf 



{T-X){T,-X) 



(9) 



We assume now that the spectrum p[X) has a spectral 
edge Ac equal to the average value of the second largest 
eigenvalue A2 of p, i.e. Ac = (A2) such that the spectrum 
for A < Ac is self-averaging. For A < Ac, close to the 
upper edge, we assume the scaling behavior 



p{X) ^ (Ac - A) 



<5s 



(10) 



that we can use to perform the integral in ([9]). More- 
over we define the spectral gap A^r of a network of size 
N as the difference between the maximal eigenvalue A 
and the spectral edge, i.e. A^v = A — Ac. Perform- 
ing a straightforward calculation under the assumption 
that the gap A^r is self-averaging in the thermodynamic 
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limit, i.e. liniTv-s-oo Ajv — ^ A, we distinguish two possible 
behaviors. If A > 0, close to the critical temperature 
T ^ Tc, we have 

^-^ = {T- T,)[l - (A^=-iC2 - Ci)/z] (11) 

where Ci_2 are constants. In this case the critical fluctua- 
tions are always mean-field. On the other hand, if A = 0, 
we have 



X 



-1 = (T - r,)[l - (T - T,)^^-'C3/z + Ci/z] (12) 



with constants C1.3. In this case the critical behavior 
depends on the particular value of 6s- For 63 > I the 
corrections of order 1/z to do not modify the crit- 
ical behavior of the susceptibility. On the contrary, for 
6s < I, the corrections of order 1/z diverge close to the 
phase transition, the fluctuations dominate the critical 
behavior and the mean-field approach cannot be applied. 
As a first check we look at the case of homogeneous de- 
gree distributions. We consider a d-dimensional lattice 
of linear size L, homogeneous hidden variables 9i = 9 \/i 
and coupling matrices 



J{n,rj) = expdf,; - fjl/do) 



(13) 



depending on the typical distance dp. In this case , we get 
always limAr_j.oo A^v = and 6s ^ [d — 2)/2 recovering 
the classical result of the Ginsburg criterion which states 
that the critical dimension for the Ising model is d = 4. 

Application to complex spatial networks - We now turn 
to the case of linking matrices p describing annealed 
scale-free networks embedded in a d-dimensional space 
with finite critical temperature Tc- For the sake of con- 
creteness, we consider a regular d = 2 lattice of side 
L and we assign to each of the N = L"^ nodes an ex- 
pected degree k according to a power-law distribution 
p{k) ^ k and we consider exponentially decaying 
couplings as in Eq. The values of the parame- 

ters {6i} are given by the solution of the set of equations 
ki = J2jPij with pij = OiOjJij. The role played by spa- 
tial embedding in the critical behavior of these networks 
is well characterized by the spectrum p{X) of the the cor- 
responding matrix p. For small do, where we expect non 
trivial effects of space, the behavior of the spectrum close 
to Ac follows Eq. ([TU)) . In Figured] we report the cumu- 
lative distribution (rank plot) of the eigenvalues of p for 
do = 1 and different values of "fsF- We observe that the 
spectral density below the spectral edge is self-averaging 
and the exponent 5s is a decreasing function of ^sf (at 
constant do) assuming values above and below = 1. 
However, the maximal eigenvalue A and the spectral gap 
A are not, in general, self-averaging, being subject to 
strong fluctuations also for large network sizes. This oc- 
curs also for the parameters values studied in Fig. [T] 
The absence of self-averaging is also observed for net- 
works without spatial embedding, where it is essentially 




FIG. 1: (Color online) Cumulative rank plot of the averaged 
spectra of 100 matrices p for scale-free random networks em- 
bedded in dimension d = 2, linear size L = 80, coupling scale 
do — 1, minimal expected degree m — 2 and jsf = 4, 6. The 
behavior for large eigenvalues is well fitted by the expression 
(|10|) with different exponent Ss below and above the critical 
value Sg = 1. Inset: average value (Aj^~^) for network en- 
sembles with the same parameters as before but with varying 
system size N — L^. The fit shows that for 73^- = 6 the 

quantity ( A^~^) increases with the system size, while for 



7SF 



4 it remains constant. 



driven by the cutoff fluctuations [21[. While this anoma- 
lous effect might be present also in spatial networks, it 
seems that the sample-to-sample fluctuations observed in 
the spectral gap are mainly due to a new feature of spa- 
tial networks, i.e. their local geometry. In fact the non 
self-averaging properties appear also for values of the ex- 
ponent 7sF (for example ^sf = 6), where the fourth mo- 
ment of the degree converges and the critical behavior as- 
sociated to complex network without spatial embedding 
is self-averaging 2l|. We checked numerically in a num- 
ber of cases that the spectral gap is non self-averaging, 
but the probability P{A^j^~^) is stable when the value 
of A^~^ is rescaled with its average value (A''^'^^) (See 
Fig. ^ . Therefore, in this case we characterize the aver- 
age critical behavior of the ensemble by the quantity 



lim 



t-t; 



lim 

W-5.0C 



A 



N 



C2-C1 



(14) 

where T^-^f is the effective critical temperature of a net- 
work and depends explicitly on the size N. If 'I' diverges, 
i.e. limjv->oo ^A^^^^ — 00, we expect that the criti- 
cal fluctuations neglected by the mean-field approach be- 
come relevant. In the inset of Figure[T]we report (A^~^) 
averaged over 100 realizations of the p matrices for the 
two network ensembles with do = 1, "/sf = 4,6 as a 
function of the network size N. The results for do = 1, 
7s_F = 4 are compatible with a limit (A^^^) — !> const 



for N ^ 00. Therefore in this case, the critical behavior 
should be well captured by the mean-field behavior. For 
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FIG. 2: (Color online) The distribution P(A'*s"^) is not self- 
averaging but is a stable distribution when the variable A*^~^ 
is normalized with its average value. Fixing the value of do — 
1, we show in plot (a) the distribution for 7s _f ~ 4, while in 
(b) 7sF = 6. Both figures are obtained from a diagonalization 
of M matrices of size N = with L = 50, 60, 70, 80. In 
particular the number of samples M is respectively M = 1000 
for L = 50, M = 400 for L = 60, M = 100 for L = 70 and 
M = 100 for L = 80. 



networks with do = I and jsf = 6, instead, (A^~^) 
seems to diverge as iV 00, signalling the presence of 
critical fluctuations not captured by the mean-field ap- 
proach. 

Conclusions - In this Letter wc have investigated how 
spatial embedding can affect the critical behavior around 
a phase transition in systems defined on spatial complex 
networks. In particular, by means of a detailed study of 
the Ising model on annealed spatial complex networks, 
we have shown that relevant critical fiuctuation not cap- 
tured by any (heterogeneous) mean-field theory may set 
in. Our analysis points out that the knowledge of the 
spectral properties of the link probability matrix p is 
crucial for the understanding of the critical behavior of 
dynamical processes and suggests a classification of the 
latter based on a generalized Ginsburg criterion. More 
precisely, when the spectrum presents a finite gap A > 
in the thermodynamic limit, the fluctuations are always 
mean-field. If instead the gap vanishes in the thermo- 
dynamic limit, the critical behavior depends on the ex- 
ponent Ss describing the scaling of the spectral density 
close to its upper edge. A fascinating open problem is 
the relation between the critical behavior of annealed and 
quenched spatial networks. The solution of this problem 
might show other new unexpected effects due to fluctu- 
ations of the local geometry. Finally our results open 
new perspectives for the comprehension of critical phe- 
nomena in spatial complex networks, whereas the general 
formalism presented here could be applied to the study 
of realistic models of epidemic spreading in transporta- 
tion networks as well as of the control of fluctuations in 
technological and power-grid networks. 
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